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EQUILIBRIUM FLUCTUATIONS FOR A DISCRETE ATLAS MODEL 


F. hernAndez, m. jara, f. valentim 


Abstract. We consider a discrete version of the Atlas model, which corresponds to 
a sequence of zero-range processes on a semi-infinite line, with a source at the origin 
and a diverging density of particles. We show that the equilibrium fluctuations of this 
model are governed by a stochastic heat equation with Neumann boundary conditions. As 
a consequence, we show that the current of particles at the origin converges to a fractional 
Brownian motion of Hurst exponent H = ^. 


1. Introduction 

The so-called Atlas model can informally described as a semi-infinite system of inde¬ 
pendent Brownian motions on M, on which the leftmost particle receives a drift towards the 
right of strength y > 0. This model is the simplest example of diffusions with rank-based 
interactions. These diffusions interacting through their rank have been proposed as a sim¬ 
ple model for capitalizations in equity markets, see (4], |7!| and the references therein. In 
0 it is proved that the equilibrium fluctuations of the Atlas model converge in a proper 
scale to the solution of a stochastic heat equation with reflection at the origin. The ultimate 
interest of this fluctuation result is that the fluctuations of the leftmost particle can be ex¬ 
pressed in the limit as a (singular) linear functional of the solution of the stochastic heat 
equation. In particular one of the main results of 0 is that the fluctuations of the leftmost 
particle are governed by a fractional Brownian motion of Hurst exponent 

A key observation in order to derive various results about the Atlas model, is that the 
sequence of spacings between the Brownian motions follows a Markovian evolution, for 
which an invariant measure of product form is known to exist. 

A natural question turns out to be what happens with a discrete analogous of the Atlas 
model. A simple discrete analogous of the Atlas model consists on a semi-infinite system 
of particles, on which the first one has a drift towards the right and subjected to the so- 
called exclusion rule: no particles can share the same position at any time t. It is well 
known that such a system of exclusion particles is in bijection with a zero-range process 
with a reservoir of particles at the origin. We learned this bijection from IfTO : although 
the system treated in CD is bi-infinite, the same bijection works for semi-infinite systems 
of exclusion particles. The earliest reference we were able to find on which this bijection 
appears is 0. In fill . Ifl3l this relation was exploited to obtain various scaling limits of 
observables of the exclusion process as a consequence of convergence results for analogous 
quantities in the zero-range process. 

For this reason, we study in this article the stationary current fluctuations of a zero- 
range process with a source at the origin. A formal description of this process is the 
following. Particles live on the semi-infinite lattice N = {1,2,...}. At each site of the 
lattice N there is a Poissonian clock of rate 2. Each time the clock at site x £ N rings, 
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one of the particles at site x moves to x — 1 or x +1 with equal probability. If the particle 
decides to move to y = 0, then the particle leaves the system. In addition, with exponential 
rate A a particle is created at site x = 1. Burke’s theorem says that for A < 1 the prod¬ 
uct measure with geometric marginals of success rate A is invariant under this evolution. 
Denote by J x (t) the signed current of particles between sites x and x+ 1 up to time t > 0. 
For x = 0, J 0 {t) denotes the number of particles that entered into the system, minus the 
number of particles that left the system up to time t > 0. The current J x (t) denotes exactly 
the displacement of the x-th particle up to time t in an exclusion process where particles 
are ordered from left to right. All but the leftmost exclusion particle are symmetric; the 
leftmost particle jumps to the right with rate 1 and it jumps to the left with rate A. In lfl3l . 
it is shown that under a diffusive scaling, the space-time limit of the fluctuations of the 
density of particles is given by a conservative stochastic heat equation, better known as an 
Ornstein-Uhlenbeck equation. With some extra work, this result can be used to derive the 
scaling limit of the current fluctuations as well. 

Let n e N be a scaling parameter. If we want the exclusion process to serve as a discrete 
approximation of the Atlas model, it is reasonable to scale A with n in such a way that the 
leftmost particle behaves like a Brownian motion with drift in the limit n -A Therefore, 
it is reasonable to take A„ = 1 — j, where b > 0 is the drift of the limiting Brownian motion. 
We will start the zero-range process with the invariant measure associated to A„, namely a 
product of geometric distributions with parameter 1 — A„. Notice that the average number 
of particles per site grows linearly with n. 

It turns out that the proof of lfl3l breaks down in that case. The heart of the proof of 
the convergence of the density fluctuation field is the so-called Boltzmann-Gibbs principle, 
which roughly states that any observable of a conservative particle system is asymptotically 
equivalent to a linear functional of the density of particles. The main issue is that the 
density of particles per site is equal to £ — 1 and it grows to infinity as n —> °°. Because of 
that, a key compactness argument in the classical proof of the so-called Boltzmann-Gibbs 
principle does not work. In j6|, a quantitative proof of the Boltzmann-Gibbs principle was 
proposed. It turns out that in our situation, this alternative proof allows to circumvent the 
compactness argument by the use of the so-called spectral gap inequality, which gives a 
sharp bound on the largest eigenvalue of the dynamics restricted to a finite box. For the 
version of the zero-range process considered in this article, the spectral gap inequality was 
proved in (l4l| . For zero-range processes with other interaction rates, see ca and na. 
According to lfl4l . the spectral gap of the zero-range process presented here has a non¬ 
trivial dependence on the total number of particles. Therefore, the proof of the quantitative 
Boltzmann-Gibbs principle of j6j needs to be adapted accordingly. 

With the Boltzmann-Gibbs as principal tool, we prove that in the stationary state the 
space-time fluctuations of the current converge to the solution of the stochastic heat equa¬ 
tion 


d r X = b 2 AX + V2W. 

with zero initial datum. In particular, as in 0 the current of particles through the origin 
converges to a Brownian motion of Hurst parameter 

This article is organized as follows. In Section[2]we introduce the zero-range process 
with a source at the origin, we define what do we understand by the current fluctuation 
field and we state the main results of this article. In Section Q] we state and give a sketch 
of proof of various estimates related to the variance of additive functionals of Markov 
processes. The exposition follows closely the one of |6), and proofs already included in 
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ED are omitted. In Section|4]we prove the main results of this article, and in Section[5]we 
discuss possible generalizations of the results proved in this article. 

2. Notation and definitions 

2.1. The model. Let £2 (l = NqO be the state space of a Markov process which we will 
describe below. We denote by 77 = { I] (x);x £ N} the elements of fin and we call them 
particle configurations. We call the elements x £ N sites and we say that 77 (x) is the 
number of particles at site x on the configuration 77 . For x,y £ N and rj £ fly such that 
?? (x) > 1 let tj xy £ £ 2 ( ) be given by 

( tj(x)- 1 ;z = x 
l ri(y) + Uz = y 
{ r)(z) -,z^x,y. 

In other words, rj x ' y is obtained from I] by moving a particle from site x to site y. For 
T) £ fly such that i](l)>lwe define 77 10 = 77 — 5] and for 77 £ £2o we define 77 0,1 = 77 + 5], 
where (>1 (x) = 0 if x ^ 1 and 5] (x) = 1 if x = 1. In other words 77 0,1 is obtained from 77 by 
adding a particle at site x = 1. 

We say that a function / : £2o —> R is local if there exists a finite set A C N such that 
whenever 77 (x) = £(x) for all x £ A, f(rj ) = /{£,). In that case we say that the support of 
/ is contained in A, and we denote it by supp(/) C A. 

Let A £ (0,1) be given. Let g : No —> R be given by g(t) = 1 if i 0 and g(0) = 0. For 
/ : flo —> R local we define Lf : flo —> R as 

Lf(r 7) = E s(v( x )) [f( r I x,x+ 1 )+f( T I x,x ~ 1 ) — 2/(77)] +A [/(T7 01 ) -/(rj)] 

for any 77 £ Qq. Although 77 (0) is not defined, it will be convenient to adopt the convention 
g(r 1 (0)) = A. Let us explain briefly Andjel’s construction of the Markov process associated 
to the linear operator L defined in this way. First we restrict ourselves to the set 

{77 g flo; E n {x)e- Mx < °°}, 

jc6N 

where M is a fixed constant. We will call this set £2 and we equip it with the product 
topology. Notice that local functions are indistinctly defined in £1 or £2 (l . We say that a 
local function / : £1 —> M is Lipschitz if there are K > 0 and ACff finite such that 

Vt 7,£ e£ 2 ,|/(Tj)— /(<u| <^E |n(*) — 5(*)|- 

xeA 

The closure of the operator L restricted to local Lipschitz functions turns out to be the 
generator of a Markov process { 77 , ;t > 0} in £1, which we call the zero-range process with 
a source at the origin. The dynamics of this process is not hard to describe. At each site 
x £ N the process waits an exponential time of rate 2 , at the end of which one particle jumps 
from x to one of this two neighbors with equal probability. If there are no particles at site x 
at the moment of jump, nothing happens. At x = 1, if the aforementioned particle decides 
to jump left, it disappears. Additionally, with exponential rate A a particle is created at site 
x = 1 . 

Let Px denote the product geometric measure on £1: 

V\(dr\) = n(l-A)A IJ Wd7j(x). 

XGN 

*We denote No = {0,1,2,... } and N = {1,2,... }. 
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In principle, Pi is defined in Oo, but it puts total mass on Cl. It can be verified that is 

_ invariant and reversible under the evolution of > 0 }. 

2.2 

2.2. The current fluctuations. Let n £ N be a scaling parameter. Let b > 0 be fixed and 
let {A n \n £ N} be defined for simplicity as A,, = 1 — j for any n > b. All the results on this 
article can be easily generalized to the case of a sequence {A n ;n £ N} in (0,1) such that 

lim n( 1 — An) = b. 

n—f °o 

Let us denote by jl" the measure p- A „ and let us denote by {r/";t > 0} the process > 

0} with initial distribution u". The time scaling may seem a little mysterious right now, 
but it will turn out to be the right one in our setting. We denote by P„ the distribution of 
0} and by E„ the expectation with respect to P„. All these parameters will be fixed 
from now on and up to the end of this article. 

For each nSN let J" (x) denote the signed number of particles crossing the bond {x,x + 
1} up to time t. Similarly, we denote by J"( 0) the number of particles created at x = 1 
minus the number of particles annihilated at x = 1, up to time t. Our aim is to obtain the 
scaling limit of the current process {/"(x);f > 0,x £ No} as n —> °°. For each function 
/ £ ^“([ 0 ,°°)) we define 

w) = E mf(-J +i E «w - p,Ml) «i) 

n «EN 0 11 xeN 

for any t > 0 and any n £ N. Here F : [0, °°) —> K is defined as F(x) = — f(y)dy for any 
x > 0 and p n = | — 1 = f 77 (x)djj." is the expected number of particles per site. 

In this way we have defined a measure-valued process {X/ 1 ;/ > 0} which we will call 
the current fluctuation field associated to the zero-range process {rj";/ > 0}. The extra 
term involving F and /] ( ',' may seem strange right now, but it will allow to get rid of a static 
drift term on the scaling limit of the current fluctuation field. 

2.3. The scaling limit. Let W be a standard space-time white noise on [0,°°) x [0,°°). We 
say that a measure-valued process {Xpt > 0} is a martingale solution of the stochastic heat 
equation 

d,X = b 2 AX + V2W (2.2) 

with reflecting boundary conditions at x = 0 if for any function / £ ^“([O.oo)) such that 
/'( 0 ) = 0 , the process 

X t (f)-b 2 f T X s (Af)ds 

Jo 

is a continuous martingale of quadratic variation 2 / f(x) 2 dx. In |3] it is proved that mar¬ 
tingale solutions of (12.2b starting from X {} = 0 are unique. Our aim is to prove the following 
result. 

| tl | Theorem 2.1. The sequence {X”;f > 0 }„<=n converges in distribution with respect to the 
uniform topology to the martingale solution of 

d r X = b 2 AX + V2W 

with initial condition Xq = 0. 

As an application of this theorem we will also prove a central limit theorem for the 
current J" (0): 


tions to a fractional Brownian motion ofFlurst index FI = 


f BM 


Theorem 2.2. The process ~hoJ" (0) converges in the sense of finite-dimensional distribu- 


def 


SHE 
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3. On additive functionals of Markov processes 

3.1. Kipnis-Varadhan inequality. The main goal of this section is to prove some esti¬ 
mates on the variance of additive functionals of the process >0}. In particular we 

want to prove the analogous of J6] Prop. 3.7], see Proposition 13.51 The proof is almost 
identical to the proof of ||6] Prop. 3.7]; we will copy the exposition of J6] Section 3.2] and 
we will explain the differences on the go. 

For f,h £ Lr{p n ) we write (ffi) = f fhdp n . we will omit the dependence in n of (•,•}, 
as well as many other quantities. Let / £ L 2 (p") such that f fdp" = 0. The //_i-norm of 
/ is defined as 

ll/ll-i = sup{2(/,/z) — (/z,—L/z)}, (3.1) 

h 

where the supremumruns over local functions in Lr{p n ). The importance of the H |-norm 
is shown by the following inequality: 


KV 


Proposition 3.1 (Kipnis-Varadhan inequality @13). For any T > 0, 


E„ 


sup 

0<t<T 




2n 



1- 


This inequality, in the form presented here was proved in J3 following the proof of a 
slightly different inequality proved in J5). 

Following E), a very efficient way to estimate the H \ -norm of a given function / can 
be achieved by using the so-called spectral gap inequality. In order to state this inequality 
we need some definitions. For £ £ N and x £ No define A^(x) = {1,... ,£} and Qp(x) = 
N^ W . For k > 0 define 


t 


= { r i e ^£^(*+0 = £}• 

i=i 


We will write Ap, £2p and instead of A^(0), Clp( 0) and 0.^(0) respectively. Let jlpj 
the uniform measure on £lp t and notice that p/, p is also the restriction of p- A to ilpj. Let us 
denote by {-,-)k,e the inner product on L 2 (pk,e)- F° r / £ L 2 (p") we denote by Var(/) the 
variance of / with respect to p". For / : Slp j -a R let Lpf : Q-p f —> R be given by 

W(t?)= E g(v(x))(f(ri x ’ y )-f(v))- 

x,y€ Ap 

We have the following proposition: 


SG 


Proposition 3.2 (Spectral gap inequality fLfll ). There exists a universal constant Kq such 
that 

(/,/) < Kq (£ + k) 2 (f-L { f) k/ 

for any k,£> 0 and any function f : Tlpp —> R such that f fdppe = 0. 


This proposition was proved in M for the zero-range process evolving on the complete 
graph and extended to finite subsets of Z d using the so called path lemma. In our one¬ 
dimensional situation, a proof can be obtained by coupling with the exclusion process. 

For x £ No and £ £ N define Ap (x) = {x + 1,... ,x + £} and 


t/(x) = 


tE^^' + O- 

1 i= 1 


varhl 
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Notice that objects like f or Lf can be defined in A fix) in a canonical way. For / : £2 —► R 
local such that supp(/) C A fix) define ysj- : Q.—> R as Yf(v) = £[/|V(jc)]. Here and in 
what follows, all conditional expectations are taken with respect to the measure fi'\ Since 
supp(/) C Af(x), the function i//| does not depend on n. We have the following proposition: 

Proposition 3.3. Let f : LI —> R be a local function such that supp(/) C A^(x) for some 
x £ No, £ £ N. Then 

II/- V'/ll-i < N/ 2 Var((1 +1] f (x))(/ — yj)). 

This proposition can be proved as Proposition 3.5 in 0. We only need to replace the 
inner products of the form (/, h) by 

*>0 

where p{k,£) = (x) = j) and make use of the spectral gap inequality on each sub¬ 

space {rj f (x) = k}. This non uniformity introduces the weighting function (1 + r/fix)) 2 
into the variance above. Taking these considerations into account, the proof of this propo¬ 
sition can be easily adapted from |6), so we omit it. 

The next proposition states that functions with supports contained on disjoint intervals 
are roughly orthogonal with respect to the H \ -norm: 


port 


Proposition 3.4. Let m £ N he given. Take a sequence 0 < xq < ■ ■ ■ < x m £ Nq and let 
{fi',i = 1,... , m) be a sequence of local functions such that supp(/) C {x,_i + 1,... ,x;} 
for any i. Define / = Xj — x,_ i. Then, 

m 

||/i+-+/m|£ 1 < KoL^Var((l+^(x ( _ 1 ))(/-v40)- 

i= 1 


Pi 


Again, the proof of this proposition is a simple adaptation of the proof of Proposition 
3.6 in 0 , so we omit it. Putting together the estimates of Propositions 13.11 and 13.41 we 
obtain the following estimate: 


Proposition 3.5. Let {f; i = 1..... } like in Proposition \3.4\ Assume in addition that 

D. 

y/f. = 0 for any i. Then 


E„ 


sup 

L o<r<r Wo “ 


m \ ' 


<^£^Var((l+T,«(fl_,))/,) 

11 i= 1 


for any T > 0. 


In what follows, Pronosition l3.5l is all we need from this section. In particular, whenever 
a estimate like the one stated in Proposition 13. 51 is available, the methods exposed in the 
following sections are independent of the spectral gap inequality, and in particular Theorem 
l2.1l holds as soon as Proposition l3.5l is available. 


intpart 


3.2. Integration by parts. For weighted differences of the functions g(rj(x)), the H i- 
norm can be estimated without appealing to the spectral gap inequality. In the context of 
hydrodynamic limits, this estimate is sometimes called the integration by parts formula. 
We have the following result: 

Proposition 3.6. Let f : No —> R be such that LvgNo /M 2 < H -00 . Then, 


xeNo 


< E /w 2 - 

xeNo 


2 
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Proof. Recall the variational formula (13.lb for the i/_i-norm. After a change of variables, 

(g(ri(x))-g(r 1 (x+l)),h) = J g{i 1 (x))(h(ri x ^ 1 )-h(n))dn n . 

Using the weighted Cauchy-Schwarz inequality and the fact that g(r)(x)) 2 = g(r)(x)) we 
see that 

2 f(x)(g(r]{x))-g{r}(x+l)),h) < Jg(q{x))(h{ri x ’ x+1 )-h{ri)) 2 d^ n + f{x) 2 . (3.2) 

Notice that 

(h-Lh)= E fg(ri(x))(h{ri^ 1 )- h (Ti)) 2 dn n . 

JteNo J 

Summing up (13.2b on x G No we obtain the desired bound. □ 


CauSch 


This estimate is not only simpler than Proposition 13.51 but also fundamental in what 
follows. The point is that the factor rj 1 (x) appearing in the general estimates has a very big 
variance compared to g(rj(x)). We will need the integration by parts formula in order to 
introduce a spatial average, after which Proposition l3.5l starts to be useful. 


Ell 


s3.1 


4. Proofs 

In this section we prove Theorems l2.ll and l2.2l The proof follows the classical structure 
for convergence theorems of stochastic processes: first we prove tightness of the sequence 
{ X" ; t > 0}„ ( 7[<j with respect to the proper topology. Then we show that any limit point of 
this sequence is a martingale solution of (12.2b with zero initial condition. Then we finish 
the proof of the convergence arguing that by the uniqueness result of 0, the sequence 
{xp-.t > has a unique limit point. 

4.1. The martingale decomposition and the continuity relation. The current processes 
Ur (x);f > 0} are Poisson compound processes with disjoint jumps. For any x G N the 
processes {M"(x);t > 0} given by 

M" (x) = J" (x) — n 4 f (g(r]"(x)) — g(T]"(x+l)))ds (4.1) 

are martingales of quadratic variation 

(M" (x)} = n 4 J‘ (gW(x))+gW(x+l)))ds. 

These formulas also hold for x = 0 if we use the convention g(J7"(0)) = A„. Since the 
currents J"( x ) have disjoint jumps, the martingales {Mf(x);f > 0}^ e N 0 are mutually or¬ 
thogonal. 

In order to simplify the notation, let us write g"(x) := g(rj"(x)). Recall the definition of 
xp(f) as a sum of currents. We have that 

K(f) =W)-W)-« 3/2 f E (g?(x)-g?(x+ l))/(f)* (4.2) 

J0 xeN 0 


martdec_vO 


martdec 


is a martingale of quadratic variation 

(K(f)) = - f E (g n s(x)+g n s(x+l))fti) 2 ds. 
n JO 


(4.3) 


variation 


Since g is bounded by 1 and E„[g”(x)] = A„, (M”(/)) is of order and it does not 
vanish in the limit n —> This observation explains the factor n 5 / 2 in the definition of 
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/: the exponent 5/2 is tied to the time scale n 4 . This still does not explain the choice of 
the time scale n 4 t. We will see that n 4 t is the time scale on which the compensator and 
the martingale part of X" (/) have the same order. We call identity (14.2b the martingale 
decomposition of Xj l (/). Using the fact that g'J(0) = A„ we can rewrite the integral term in 
(14.2b as 

rt 


x<EN 

where V"/ := n(/(~) — is a discrete approximation of /'(A). 

Apart fromx = 1, particles are not either created nor destroyed by the dynamics. There¬ 
fore, for x > 2 we have the relation 


» 1/2 / E {s n s(x)-K)Kfds, 


(4.4) 


int 




(4.5) 


cont 


In other words, what comes in minus what comes out is equal to what we have minus what 
we had. For x = 1 this relation, properly understood, also holds. We call identity (14.5b 
the continuity relation. This relation will allow us to express the integral term in (14.2b in 
terms of the process X ", aside from an error term that goes to 0 as n —> Let us explain 
in a heuristic way how will we do this. Notice that g”(x) is a Bernoulli random variable 
and f (1 — g(ri(x)))dp n = 1 — A„ is asymptotically equivalenQto 4-, where p n = | — 1 = 
f r\(x)djl". The so-called Boltzmann-Gibbsprinciple states that at the level of fluctuations, 
any local function can be approximated by a function of the density of particles. By the law 
of large numbers, this function should be A in first approximation. Since we are looking 
at these variables at the level of the central limit theorem, it is reasonable to assume that 
the average of g(p(x)) — A„ on the box Ai(x) is well approximated by 1 ■, (p 1 (x) — p„). 

j_J~Pn 

If we replace g"(x) — A„ by . 1 (tj"(x) — p n ) in the integral term of (14.2b and we use the 
continuity relation (14.5b . we can rewrite this integral term as 

b 2 r E J n s WA/(|)^ = b 2 [x?(Af)ds 


plus error terms. It is exactly at this passage that we need to use the hypothesis /'(0) = 0. 
Otherwise a non-vanishing boundary term would appear. Therefore, assuming the validity 
of the Boltzmann-Gibbs principle it is possible to rewrite ( 14.2b as an approximate closed 
equation for the current fluctuation field X,". 

Notice that the function g"(x) is most of the time equal to 1, being only eventually 
equal to 0. Therefore, the dynamics is very singular and the Boltzmann-Gibbs principle 
can not be proved by the traditional method introduced in m (see iflOl Chapter 11] for 
a comprehensive reference). For this reason we need to use the quantitative proof of the 
Boltzmann-Gibbs principle introduced in |6 |. 

4.2. Tightness. In this section we prove tightness of the sequence of processes {A,"; t > 
0}„ 6 fj. As usual, we restrict ourselves to a finite time horizon [0,7’]. Tightness for the 
process in [0,°°) follows pasting intervals of fixed size T. 

Recall that we are thinking about {X' t u ,t € [0,7’]} as a measure-valued process. This is 
very convenient, since we can reduce tightness considerations to real-valued processes: 


^ Loosely speaking, we will say that two sequences a n , b„ are asymptotically equivalent if j 2 - —> 1 as n —> 

To be precise, when the sequences are composed of random variables, we should specify the sense on which the 
limit holds, but this point will not be of any relevance. 
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Proposition 4.1. The family { X "; t £ [0.7'] }„cp; is tight with respect to the uniform topol¬ 
ogy if and only if for each function f £ the family of real-valued processes 

{X"(f);t £ [0,7’]}„ eN is tight. 


martwhitt 


This proposition is standard in the context of hydrodynamic limits of interacting particle 
systems. A proof of it on finite volume can be found in [TO] Chapter 4], The proof adapts 
easily to the unbounded case. 

Recall the martingale decomposition (14.2b . The proof of tightness for X"(fj can be 
reduced to the proof of tightness of the martingale M' r \f ), the initial distribution Xq (/) 
and the integral term in (14.2b . A simple computation shows that the initial distribution 
Xq (/) converges in distribution to a Gaussian random variable of mean 0 and variance 
f F(x) 2 dx. This integral is finite since F £ ^“([0,°°) by definition. 

For martingales, powerful methods are available. The integral term will be more de¬ 
manding, but Kipnis-Varadhan inequality (13.1b coupled with the integration by parts stated 
in l3.6l will provide the necessary bounds. 

Let us state a convergence criterion for martingales, see □3 Theorem 2.1] 

Proposition 4.2. Let {Mf;f £ [0,r]}„ e fj be a sequence of martingales with M'q = 0 and 
let A'j. the size of the biggest jump ofM'j in the interval [0,7’]. Assume that 

i) (M”) converges in distribution to <J 2 t, 

ii) A" converges in probability to 0. 

Then {M'/\t £ [0,7’]}„ 6 pj converges in distribution to a Brownian motion of variance G 2 . 


The integral term will be handled with Kolmogorov-Centsov’s tightness criterion (see 
J8] Exercise 2.4.11]): 


KC 


Proposition 4.3 (Kolmogorov-Centsov’scriterion). Let {Yf\t £ [0.7 ’]} n cf \ be a sequence 
of real-valued processes with continuous paths. Assume that there exist constants K,a,a' > 
0 such that 

E[\Y t n -Yf\ a \<K\t-s\ l+a ' 

for any s,t £ [0,7’] and any n £ N. Then the sequence { Yf ; t > 0} ne N is tight with respect 
to the uniform topology. 


Now we are in position to prove the tightness of {X";t > 0}. Notice that the current 
processes Jj (x) have jumps of size 1. Therefore the jumps of Mj (fj are at most of size 
In particular the martingales Mj (fj satisfy part ii) of the convergence criterion. 
Recall the martingale decomposition (14.2b and the formula (14.3b for the quadratic variation 
of M”(/). In order to prove i), it is enough to observe that 

fimE n [{M\\f))] =2t\\f\\ 2 L2{R) 

and that 

[( (K (/))- KK (/))]) 2 ] < % E/(f) 2 

n *eN 

for some constant C depending only on {A„}„ g n. Therefore, not only the martingale se¬ 
quence {M,"(/);/ £ [0,7’] is tight but it also converges to a Brownian motion of vari¬ 
ance 2 f f(x) 2 dx. 

The integral term ( 14.4b is more demanding. Let us introduce the definitions 

* W 0 






10 


F. hernAndez, m. jara, f. valentim 


_ g n s{x+\) + --- + g n s {x+t) 

8s W £ 

Let /; : No —> R be such that Jl x h(x) 2 < +°°. Notice that || • ||_i satisfies the triangle in¬ 
equality. Using the triangle inequality twice we see that 


E (*(»/(*))-^WM *)||_1 <71111 E {8( r l(x+i-l))-g('n(x+i)))h(x)\ 


xGN 0 


;=l<=l N 0 


< , 


( E h ( x Y 

jteNo 


1/2 


Combining this estimate with Proposition l3.1l we obtain the bound 

>2-1 mf 


E„ 


(» 1/2 / E (*;(*)-^(*)) v s/<fo) 

J0 agNq 


< 


E ( v "/) 2 > 


(4.6) 


AGNo 


which is of order By Cauchy-Schwarz inequality. 


E„ 


(» 1/2 f E (g^w-^Kfds) 2 } < bJ j- E ( v */) 2 

J0 aGNq 1 AGNq 


which is of order Choosing i = [nf *( 3 ] we have just proved that there exists a 

constant C : = C(f) such that 


E „ 


U' 2 f E («?(*) 2 1 <ct 5/3 

V J0 agN ' J 


for any n € N and any t g [0, T], Since the increments of this process are stationary, we 
have just proved that the hypothesis of Pronosition l4.3l holds for the integral term (14.4b with 
a = 2 and a' = |. In consequence, the processes 

« 1/2 / E (<?" M - K) Kf ds 

JO r<=M 


are tight. We conclude that {X"(f)\t G [0,r]}„ G ^ is tight for any / G ^TdP, 00 )) an d by 
Proposition ^. 1 I the measure-valued processes {X";t g [0, are tight. 


4.3. The Boltzmann-Gibbs principle. Our objective in this section is to prove the so- 
called Boltzmann-Gibbsprinciple for the integral (14.41 ): 


BG 


Proposition 4.4 (Boltzmann-Gibbs Principle). For any function f g 'i£’, oo ([0. <*>)), 


lim E„ 


,/u jcgN 


2-i 


= 0 . 


(4.7) 


What this proposition is telling us is that the integral (14.41) is well approximated as n —> °° 
by a linear function of the density of particles. The proof of this theorem is somehow 
winding. We will successively prove that the integral (14.41 ) is asymptotically equivalent to 
other expressions as n —> until we end up with the density of particles. First we will 
introduce a spatial average on g"(x). Then the main step comes, which is to replace spatial 
averages of g" (x) by a function of the particle density 


V?(x) 


v"(x+l) + --- + V ” 0 + 1 ) 

£ 


estl 


BGliml 
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Then we show that this function of rj' Ll (x) is well approximated by its linearization around 
p„. Finally we undo the spatial average to recover the required estimate. 

The following lemma is just a slight modification of estimate (14.61) . so we state it without 
proof. 


| 11 | Lemma 4.5. For any function f £ c tf’f‘([0,°°)), 

K\U /2 f E {g n s (x)-g n /{x-\))V' x fds\ 
jo TgN 


< ^ E (v?/) 2 . 

xeNo 


What this lemma has accomplished is to replace g n s (x) — X n by g n s ' ( (x— 1) — A„; the latter 
has a spatial average of size l. Notice that the sum g"’ (x— 1) starts at x. For l £ N define 
y/ : Q —» R as 

V x (n)=E[g(Tl(x+i))\ri e {x)]. 

An explicit computation shows that 


For rSN and t > 0, define 


Vrf(rj) = l- 


V?’ e (x) = &&?). 


The core of the proof is the following lemma: 


(4.8) 



Lemma 4.6. For any function f £ «T([o, “)), 



where C(f) is a constant which depends only on f and the parameters of the model. 


Proof Notice that if \x — x’\ > £ then the supports of the functions 

^(x)-^(x), g^(x')- Y ^(x') 

are disjoint. Therefore, at the price of a multiplicative constant ( we can put ourselves into 
the setting of Proposition [33] Therefore, the expectation above is bounded by 

18 y E ( V "+i/) 2 Var (l 1 + V e (*)) (/« - V(*)))■ (4-9) 

n j:6No 

In order to estimate the variance in ( 14.91 ) we use the elementary inequality 
Var(AT) < [e(X- P x ) A ] 1/2 [EyY 2 + 2px[EY 4 ]^ 2 [E(X - p x ) 2 ] 1/2 + p^EY 2 , (4.10) 
which comes from the identity 

X 2 Y 2 =(X- p x ) 2 Y 2 + 2p x Y 2 (X - p x ) + p 2 Y 2 , 

valid for any random variables X, Y with means px and zero, respectively. 

Taking X = 1 + rj^(x) and Y = g e (x) — y e (x) we see that EX = E(X — p x ) 2 < 
and EY 2 < Furthermore, there exists a finite constant C such that E(X — p x ) < ( -jr 
and EY 4 < Putting this estimates into (14. 101) . we obtain the bound 

Var((l + (x)) (g\x) — y/(x))) < 


explicit 


boundl 


var 

















12 


F. hernAndez, m. jara, f. valentim 


and we conclude that (14.9b is bounded by 


Ctl 

n 3 / 2 


I(v;/) 2 , 

xeN 


which completes the proof. 


□ 


Proof of Proposition 1 4.4\ Let us summarize what we have done up to here. On one hand, 
combining Lemmas 14.51 and l4~6l we see that the integral term (14.4b is asymptotically equiv¬ 
alent to 

n 1/2 f E ('K'\x)-K)V n x+1 fds 

jtGPj 

as soon as £ -C h 1 / 2 . On the other hand, it is easy to check that the term T]"(x) appearing in 
the formulation of the Boltzmann-Gibbs principle can be replaced by the spatial average 
rj"' e (x) whenever £ <C n. Therefore, in order to prove Pronosition l4.4l it remains to verify 


lim E„ 


(» 1/2 f E (K’'W-4-( I ^ I (ri^(x)-p n ))V” +I /^) 

■'0 xGN 


2n 


= 0. (4.11) 


The same argument used at the beginning of the proof of Lemma [47)1 permits to bound 
the preceding expectation by 

. 2"i 


Cfn 2 / 2 E„ 




(4.12) 


where C = C(f) is a constant depending solely on f. 

Observe that the expression squared into the expectation above almost corresponds to 
the error committed in the linearization of the function h(z) = (1 +z ) -1 around p„ evalu¬ 
ated at T) n ’ e (x), which is explicitly given by 

1 1 (z~P») = (■ z-p n ) 2 

1+Z 1+Pn (1+Pn) 2 (1 +z)(1+Ph) 2 


The word almost in the preceding paragraph is due to the term j appearing in d4~8b . 
Rearranging terms in a convenient way, the expectation in (14.12b can be bounded above by 
twice 


1 

(1 +P«) 4 


E„ 


' (l] n,e (x) — Pn) 4 ' 
. {l + ri"’ e (x)) 2 . 


1 

ww 


■E„[(l + J7 '%))- 2 ]. 


(4.13) 


Given a > b, let us define a„ = | — 1 Considering separately the cases {t]”’ ( (x) > a„} and 
{77"’ f '(x) < fl„} we can bound the two expectations above by 


4e„ 


(tl"’ e (x)-p„) 4 +p£ n"(TJ n ’ e (x)<a n ) and 4 +fl n (T] n ’ e (x) <a„), (4.14) 


respectively. 

From ( 14.12b -( 14.14b and the fact that ¥. n \(r) e (x)—p n ) 4 ] = / ft), we see that in order 

to prove (14.11b it is enough to show that 

a 2 (j? + jjzjyi) +^n 2 { 1 + {i \ )i )p n ( T ] n ’ e {x) < a n ) (4.15) 

goes to zero as n goes to infinity. According to (14.16b below, the probability p n (rj n ’- (x) > 
a„ ) decays exponentially fast in £. Therefore, the expression in (14.15b goes to 0 as soon as 
1<L This ends the proof of the Boltzmann-Gibbs principle. □ 


BGlim2 


eh 


e2 


e3 


ec4.1 
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current 


For the sake of completeness, we conclude this subsection with the derivation of the 
exponential estimates mentioned above for tail probabilities of Tj e (x). The estimates are 
not completely standard due to the increasing density of particles. 

Let X be a random variable with distribution Geom(0): geometric distribution of suc¬ 
cess probability 6. Notice that p := E[X\ = Furthermore, denoting by .#p(A) := 
E[e Xx ] the moment generating function of X, we have that 


„#p(A) 


1 

1 -p(e A -l) 


for A < ln(ii£) and .# p (A) = +°° otherwise. 

Notice that r\ n '' i (x) is the average of l independent geometric random variables with 
success probability |, thus p„ := E[q n (x)\ = | — 1. Therefore, Cramer’s method allows us 
to obtain the following exponential bounds on tail probabilities: for any a > 0 

\\ogP(r\ n/ (x) <a)< (4.16) fcT| 


where . Ap denotes the large deviations rate function associated to geometric distributions 
of mean p: 


J ! p(a) := sup {Aa —log^p(A)} 

AeR 


a log 


a(l + p) 
p (1 + a) 


log 


1 -f- Cl 

1 +p' 


On the other hand, it is no difficult to see that taking a n = | — 1 we have 


lim J^p («n) = 5 — log jj — 1 . 

n —ino ' M *'* 


(4.17) c2 


Observe that the right hand side of the last line coincides with the large deviations rate 
function associated to exponential distributions. Indeed, this is consistent with the well 
known fact that if X" has distribution Geom(^), then jX n converges to an exponential 
distribution of mean jj. 


4.4. The convergence. Now we are in place to prove Theorem 12.11 In Section 14.21 we 
showed tightness of the sequence {X";t £ [0, T] }„ 6 n with respect to the uniform topology. 
Then there are subsequence n! and measure-valued process {X™\t £ [0,7]} with continu¬ 
ous paths such that {X"';t £ [0,7"]} converges to {X t °°;t £ [0,7]} in distribution with respect 
to the uniform topology. Recall the continuity relation stated in (14.5I >. We have the relation 


W 3 / 2 


E 


(r?;W-pn)V"/ 


4n E E (Vo(*)-Pn)Kf 

n / xeN 0 n / xm 


(4.18) 


cont2 


where A"/ := u(V" +| / — V"/) is a discrete approximation of A/(£). Notice that the last 
term on the right-hand side of this identity is a boundary term. The following proposition 
will prove to be useful. 


Proposition 4.7. Let h : No —>• R be such that £ a . gNq h(x) 2 < +°°- Then, 

, 2-i 


E„ 


^ 52 < 32 tn 4 52 h{x) 2 . 


jcGNq 


*6No 
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Proof. Multiplying j4.lt by f(x) and adding up on x we obtain the martingale decomposi¬ 
tion 

E J ' r \x)h(x) = E M?(x)h{x)+n 4 ( E {s n s W) -g n s {x+\))h{x)ds. 
agNq aGNq agNq 


Recall that the martingales M"(x) are mutually orthogonal. Therefore, 

2n 


E„ ^ E Af" (x)h(x)j < 2n 4 t E h(x 


AG No agN 0 

The integral term can be estimated by combining Propositions 13 .1 1 and 13.61 We obtain the 
bound 


E„ 


(« 4 / E (8 n s(x)-g"(x+l))h{x)ds) 

~ v JQ agN 0 

By using the triangle inequality, 

£[(W + 7) 2 ] < (E[X 1 } l l 2 + E[Y 2 ] 1 l 2 ) 1 


< 18fn 4 E M x ) 

agNq 


the proposition is proved. 


□ 


This lemma is useful for two things. First we see that 
b 2 




agN 0 


2 1 Vch 4 

<--EM))' 


H 


agNq 


and in particular the first sum on the right-hand side of equation (14. 1 8b is equal to 

t?XXAf)-*Z t (T t S{x)-p n )f{z) 


aGN 


plus an error term that vanishes in L 2 (P„). The second sum on the right-hand side of (14. 1 8b 
is equal to 

4t2 E (VoM-Pn)/(l) 

n agN 

plus an error term that vanishes in L 2 (P n ). Therefore, the two first sums on the right-hand 
side of (14.18l > are equal to b 2 Xp(Af) plus an error term that vanishes in L 2 (P„). It is exactly 
here that we need the extra sum on the definition of xp(f) introduced in (12. 11 1 . 

The second use for Pronosition l4.7l is to show that the third term on the right-hand side 
of (14. 1 8b is negligible. We have that 

E„[(n- 3 / 2 T s "(0)V”/) 2 ] < 32sn(Vg/) 2 . 

It is exactly at this point that we need to assume that /'(0) = 0. In that case, |Vq/| < 
^||A/||oo and the boundary term in (14. 1 8b vanishes in L 2 (P„) as n —> By the Boltzmann- 
Gibbs principle stated in Proposition 14.41 we conclude that 

« 1/2 / E {a1{x)-Ki)y n x fds = f Xp(Af)ds 

JO Jo 


plus an error term that vanishes in L 2 (P„) as n —> In particular, (14.21 ) can be written as 

K(f) = W) -Xo (/) - f K{b 2 Af)ds 

Jo 
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plus a term that vanishes in Zr(P„). Taking limits along the subsequence n' we conclude 
that for any / 6 ^“([0,°°)) such that /'(0) = 0, 

*T(/)-W)- fx?(b 2 Af)ds 
Jo 

is a Brownian motion of variance 2 / f(x) 2 dx. Recall that Xfi (/) converges to a Gaussian 
random variable of mean 0 and variance p / f{x) 2 dx. Therefore X^(f) is a spatial white 
noise of variance Jy. In other words, {A - ”;/ £ [0, T’]} is a stationary solution of the sto¬ 
chastic heat equation (12.2I > and it is unique in distribution. We conclude that the sequence 
{X";t £ [0, r]}„ e pj has a unique limit point and therefore it converges to it, which proves 
Theorem l2.ll 


4.5. Proof of Theorem 12.21 The proof of Theorem [Z2] follows an idea which we learned 
from |16l . Let (p : (0,°°) —> [0,°°) be a smooth function of support contained on (0,1), 
such that f (p(x)dx = 1. For e > 0 define (p £ (x) = |<p(|) and h e (x) = ^(p(^)dy. Notice 

that h £ (0) = 1, h £ (x ) = 0 if x > £ and {(p £ ;e > 0} is an approximation of the 5 of Dirac. 
Using the continuity relation ( 14.5b we see that 

-4 too) - ETON© -*.(^») = -4 f WM - 

n x x=l ' n x=l 


Up to some small error term, the sum on the left-hand side is just Xj\(p £ ). In particular, 
there is a constant C depending only on the parameters of the model and the choice of <p 
such that 


E n 


(i«°)-*"(«><)) 


2i 


Taking n —> °° and then e 


L 

0 we conclude that 


< Ce. 


lim 

oo 


172^(0) = limX( 9e ) 


in the sense of finite-dimensional distributions, as soon as the right-hand side is well de¬ 
fined. As observed in 0, this limit exists and it is equal to a fractional Brownian motion 
of Hurst exponent H = which proves Theorem 12.21 


5. Discussion and generalizations 

Non nearest-neighbor transition rates. As mentioned in the introduction, there is a com¬ 
binatorial relation between the zero-range process studied in this article and the exclusion 
process with symmetric particles with the exception of the leftmost one. The large-density 
limit considered in this article corresponds to a vanishing density of particles in the exclu¬ 
sion process. Therefore, a proof of the main result of this article using the exclusion process 
representation and following the steps of 0 does not seem to be out of reach. Neverthe¬ 
less, the proof presented here has the advantage of being more general, since it is built upon 
general properties of interacting particle systems, namely the spectral gap inequality, the 
equivalent of ensembles and the product structure of the invariant measure. Let us mention 
here a simple generalization of the model for which our proof can be adapted. Let p( ■) be 
a symmetric transition rate in Z with finite range, but not necessarily equal to 1. For y < 0, 
x > 0 and q e Q such that I] (x) > 1 we define r) x - y = rj - 8 X . The operator 

Lf(V) = E p(y - x)lMx) > 1} [/(q^) - /(q)] + A E p( x ) [Av + « - /(r?)] 

XGN XGN 

y£ Z 
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defines a zero-range process in Q which has the same invariant measure as the zero-range 
process with nearest-neighbor jumps introduced in Section [2] Let us use the same notation 
> 0} for the process generated by L. Currents need to be defined in a different way. 
Let us denote by R the range of p(-). For x,y £ N with |y — x\ < R we define J'/(x,y) the 
signed number of particles passing from x to y and let us define 7"(0.x) as the number of 
particles created at x minus the number of particles destroyed at x. Let us define 

E Jt^y). 

0<z<x 

y>x 

The current fluctuation field X" is defined like in Section 12.21 using this version of the 
current processes. Then Theorem lXTl holds as well: 

Theorem 5.1. The sequence { X "; t > 0} ng pj converges in distribution with respect to the 
uniform topology to the martingale solution of 

d,X = b 2 0 2 XX + V2# 

with initial condition Xq = 0, where <J 2 = Lr>0 <- 2 p(z). 

The interested reader will not find any difficulty to adapt the proof of Theorem 12. II for 
this case. The only relevant difference is the following version of the continuity equation 
(l431 >: for any x < y such that y — x>R, 

t W(z)-ril ?(*))• 

z=x+ 1 

General interaction rates. As mentioned before, the main properties of the zero-range 
process used in the proof of Theorem |5T| are a sharp bound for the largest eigenvalue of 
the dynamics restricted to a Ante box, the explicit knowledge of the invariant measure of 
the system and the equivalence of ensembles. With these three properties in hand, the 
proof follows through. In fl5l . the spectral gap inequality for the zero-range process with 
interaction rate g(k) = kf, y (0,1) was derived, as well as the estimates needed to obtain 
the equivalence of ensembles. Therefore, the proof of Theorem[2T]could in principle be 
adapted to the zero-range process with interaction rate g(k) = kf. The limiting equation 
will be the same of Theorem 12. II with diffusion coefficient and noise variance depending 
on the parameters of the model. 

Acklowledgements 

F.H. thanks FAPERJ for its support through the grant E-26/112.076/2013. M.J. thanks 
CNPq for its support through the grant 401628/2012-4 and FAPERJ for its support through 
the grant JCNE E17/2012. M.J. would like to thank Amir Dembo for valuable discussions. 
F. V. would like to thank the warm hospitality of IMPA, where this work was done. 

References 

[1] Chih Chung Chang. Equilibrium fluctuations of gradient reversible particle systems. Probab. Theory Related 
Fields , 100(3):269-283, 1994. 

[2] Chih-Chung Chang, Claudio Landim, and Stefano Olla. Equilibrium fluctuations of asymmetric simple 
exclusion processes in dimension d > 3. Probab. Theory Related Fields, 119(3):381—409, 2001. 

[3] Amir Dembo and Li-Cheng Tsai. Equilibrium Fluctuation of the Atlas Model. ArXiv e-prints (2015). 
arXiv:1503.03581 

[4] E. Robert Femholz. Stochastic portfolio theory, volume 48 of Applications of Mathematics (New York). 
Springer-Verlag, New York, 2002. Stochastic Modelling and Applied Probability. 


EQUILIBRIUM FLUCTUATIONS FOR A DISCRETE ATLAS MODEL 


17 


[5] Pablo A. Ferrari. The simple exclusion process as seen from a tagged particle. Ann. Probab., 14(4): 1277— 
1290, 1986. 

[6] Patricia Gonsalves and Milton Jara. Scaling limits of additive functionals of interacting particle systems. 
Comm. Pure Appl. Math., 66(5):649-677, 2013. 

[7] Tomoyuki Ichiba, Vassilios Papathanakos, Adrian Banner, Ioannis Karatzas, and Robert Femholz. Hybrid 
atlas models. Ann. Appl. Probab., 21(2):609-644, 2011. 

[8] Ioannis Karatzas and Steven E. Shreve. Brownian motion and stochastic calculus, volume 113 of Graduate 
Texts in Mathematics. Springer-Verlag, New York, second edition, 1991. 

[9] C. Kipnis and S. R. S. Varadhan. Central limit theorem for additive functionals of reversible Markov pro¬ 
cesses and applications to simple exclusions. Comm. Math. Phys., 104(1): 1—19, 1986. 

[10] Claude Kipnis and Claudio Landim. Scaling limits of interacting particle systems, volume 320 of 
Grundlehren der Mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences]. 
Springer-Verlag, Berlin, 1999. 

[11] C. Landim, S. Olla, and S. B. Volchan. Driven tracer particle in one-dimensional symmetric simple exclu¬ 
sion. Comm. Math. Phys., 192(2):287-307, 1998. 

[12] C. Landim, S. Sethuraman, and S. Varadhan. Spectral gap for zero-range dynamics. Ann. Probab., 
24(4): 1871-1902, 1996. 

[13] Claudio Landim and Sergio B. Volchan. Equilibrium fluctuations for a driven tracer particle dynamics. 
Stochastic Process. Appl., 85(1): 139—158, 2000. 

[14] Ben Morris. Spectral gap for the zero range process with constant rate. Ann. Probab., 34(5): 1645-1664, 
2006. 

[15] Yukio Nagahata. Spectral gap for zero-range processes with jump rate g(x) — xf. Stochastic Process. Appl., 
120(6):949-958, 2010. 

[16] Hermann Rost and Maria Eulalia Vares. Hydrodynamics of a one-dimensional nearest neighbor model. In 
Particle systems, random media and large deviations (Brunswick, Maine, 1984), volume 41 of Contemp. 
Math., pages 329-342. Amer. Math. Soc., Providence, RI, 1985. 

[17] Ward Whitt. Proofs of the martingale FCLT. Probab. Surv., 4:268-302, 2007. 

INSTITUTO DE MATEMATICA, UNIVERSIDADE FEDERAL FLUMINENSE, RUA MARIO SANTOS BRAGA 
S/N, Niteroi, RJ 24020-140, Brazil 
E-MAIL: freddyhernandez@id. uff. br 

IMPA, Estrada Dona Castorina 110, Rio de Janeiro, RJ, 22460-320, Brazil. 

E-MAIL: mjara@impa.br 

DEPARTAMENTO DE MATEMATICA, UNIVERSIDADE FEDERAL DO ESPIRITO SANTO, AV. FERNANDO 

Ferrari, 514, Goiabeiras, Vitoria, 29075-910, Brazil. 

E-MAIL f abio. valent im@ufes . br 



